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ON THE PARABOLICITY OF THE MUSKAT PROBLEM: 
WELL-POSEDNESS, FINGERING, AND STABILITY RESULTS 

JOACHIM ESCHER AND BOGDAN-VASILE MATIOC 



Abstract. We consider in this paper the Muskat problem in a periodic geometry 

and incorporate capiUary as well as gravity effects in the modelling. The problem 

rr\ ' is rewritten as an abstract evolution equation and we use this property to prove 

well-posedness of the problem and to establish exponential stability of some flat 
equilibrium. Using bifurcation theory we also find finger shaped steady-states which 
are all unstable. 



1. Introduction 



< 

C^ , The Muskat problem was introduced into the mathematical literature by M. Muskat 

[19] as a model for the simultaneous flow of two immiscible fluids in a porous medium. 
The problem is of considerable practical relevance since it is also a model for the 
interaction between oil and water in an oil reservoir. In practice high flow water 
^ ' injection pumps fill the bottom of the reservoir with water to push the oil towards the 

CN , wells like a piston. Furthermore the injected water is used to keep the pressure in the 

. —^ I reservoir unchanged over a long term, which is keeping the production rate constant. 

^v^ ■ From the practical point of view it may be of interest to determine the optimal speed 



at which water drives the oil, which adheres more strongly to the porous structure, 

Q ■ upwards to the oil well. 

^D . In contrast to the Mullins-Sekerka problem, where the interface separating the two 

fluids is driven forward with normal velocity equal to the normal component of the 
velocity jump over the interface, cf. [12], the Muskat problem considers a continuous 

K> , velocity field in the whole domain occupied by the two fiuids. The authors in [12] show 

^ ' that the Mullins-Sekerka problem has a very natural parabolic structure and solve it 

by using the rich literature on this topic. A first result on the classical solvability 
of the Muskat problem was obtain in [22], where local in time existence of classical 
solutions is established by using Newton's iteration method. Later on, in [23], the 
Schauder fixed point theorem and existence results for a elliptic diffraction problem 
are used to prove existence of global solutions for small perturbations of some fiat 
interface. 

The problem has been also considered by means of complex methods. We like to 
mention the work [21] where global in time existence for initial data that are small 
perturbations of a fiat interface are proved. The authors of [21] also show that the 
problem is ill-posed when the less-viscous fiuid drives forward the fiuid with higher 
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viscosity. Further local well-posedness results have been established [1] and more 
recently in [5,6] by means of energy estimates in the absence of surface tension forces. 

Regarding the stability of equilibria an interesting result is found in [13] where the 
stability of some circular steady-state is proven by considering series expansions and 
surface tension effects at the boundary between the fluids. The authors of [13] state 
that the equilibrium is not generally asymptotically stable. As far as we know there 
are no other stability results for the Muskat problem. Both, the Mullins-Sekerka and 
the Muskat problem are also closely related to fingering phenomena which have and 
still receive a lot of interest in many fields of sciences. 

On the parabolicity of the Muskat problem not many results are available. We show 
in this paper that the Muskat problem and the Mullins-Sekerka model are related 
in the sense that the Muskat problem is, at least in a neighbourhood of some flat 
interface, of parabolic type too. This observation is true, when considering surface 
tension effects, independently of the boundary data. When allowing only for gravity 
and viscosity effects we need to make restrictions to the boundary data which are 
consistent with earlier results mentioned above. The approach proposed here is not 
only a useful tool to prove the well-posedness of the problem but enables us also to use 
the theory on parabolic problems to establish the dependence of the solutions on the 
initial data (see Theorem 2.1). Furthermore, parabolic theory provides us a simple 
argument to prove exponential stability of some flat equilibrium, cf. Theorem 5.3. 
Additionally, Corollary 2.4 yields an optimal value for the normal velocity at which 
water may displace oil in a oil reservoir. When considering surface tension effects, 
we find in Theorem 6.1 for small values of the surface tension coefficient steady-state 
fingering solutions which are all unstable, cf. Theorem 6.3. In the Appendix we show 
that the problem under certain constant boundary conditions corresponds to a Muskat 
problem in a frame moving with constant velocity. 

The outline of the paper is as follows. We start with a short description of the model 
of interest and present the well-posedness result in Section 2. In order to establish 
this result, we transform in Section 3 the problem on a fixed domain by straightening 
the unknown boundary and reduce the problem to an abstract evolution equation. In 
Section 5 we prove the asymptotic stability of some fiat equilibrium which appears 
under certain constant boundary conditions. In the following section we find, using 
the surface tension coefficient as a bifurcation parameter, global bifurcation branches 
consisting only of unstable steady-state solutions of the problem. In the Appendix we 
analyse the Muskat problem in a moving frame. 

2. Mathematical model and the well-posedness result 

To give a precise description of the physical setting we are interested in, consider 
a porous medium (or a vertical Hele-Shaw cell) containing two incompressible and 
immiscible Newtonian fluids in motion. We use the subscript — when referring to the 
fluid located in the lower part of the porous medium and which occupies the domain 
ri_(t), respectively -|- for the other one located in r2+(i). These two fluid phases are 
separated by the interface r(t) which evolves in time and is to be determined as a part 
of the problem. The physically relevant problem is essentially three-dimensional, but 
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n_(t) 



Figure 1 . The 2-D setting 



it may be approximated by a two-dimensional mathematical model. We refer to [15] 
for a deduction of a generalised Darcy's law for non-Newtonian fluids. The methods 
presented there apply also to Newtonian fluids and yield for a fluid in motion in a Hele- 
Shaw cell the well-known linear Darcy's law for a two-dimensional averaged problem. 
This is the reason why we restrict our considerations to the situation ^±{t) C M^. 
The Muskat problem is a potential flow in the sense that the velocity fields v± satisfy 
Darcy's law 

k 

v± = \/u± in ri-i-(i), (2.1a) 

/^± 
cf. [19]. The potentials 

u±:=p±+gp±y in n±{t) (2.1b) 

incorporate both the dynamic pressure p± inside i^±{t) as well the hydrostatic term 
gp±y- We write p± to denote the densities of the fluids, g is the gravity constant, k 
the permeability of the porous medium, p± are the viscosity constants, and we have 
chosen y to be the height coordinate (see Figure 1.). Incompressibility means that 

\J.v±=Q in VL±{t). (2.1c) 

The pressure jump across the interface T{t) is compensated by the surface tension in 
the interface. From the Laplace- Young condition we then obtain 

u+-u- = 7Kr(i) + g{p+ - p-)y on r(t), (2. Id) 

where 7 is the surface tension coefficient of the interface, and k-yU) is the curvature 
of r(t). Two more, so-called kinematic boundary conditions on T{t) are given (see 
also [13,19,22,23]) by 

V{t,-) = {v±{t,-),v{t,-)) on r(t), (2.1e) 

meaning that the interface moves along with the fiuids (a particle on the boundary 
remains on the boundary as time elapses). Here, v{t) stands for the unit normal at T[t) 
which points into 17+ (t). In order to study the problem we are left to impose boundary 
conditions on the fixed interfaces r_i and Fi. On the bottom r_i we presuppose the 
value of U- to be known, 

u- = gi{t,x) on r_i, (2. If) 
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and on Fi we set 

df7U+ = g2(t,x) on Ti, (2.1g) 

where V = (0, 1) is the unit normal vector at Ti which points outward from i7_|_(t). Of 
course, it is also possible to choose Dirichlet conditions on r_|_ and Neumann conditions 
on r_i, but we prefer to remain within the setting considered in [22,23]. The interface 
at time i = is also prescribed: 

r(0) = To. (2.1h) 

System (2.1) is a two-phase moving boundary problem. The main interest is in 
determining the motion of the interface r(t) separating the fluids. If T{t) is known, 
then the potentials u± can be then found by solving elliptic mixed boundary value 
problems. We shall prove that if the boundary Fq is the graph of a function which is 
small in some suitable Banach space, then, locally in time, problem (2.1) possesses a 
unique classical Holder solution. To that aim, we first rewrite problem (2.1) in a more 
accessible way by introducing a parametrisation of the, a priori unknown, moving 
boundary T(t). 

2.1. Parametrising the boundary. For simplicity, we consider in the following horizon- 
tally 27r— periodic flows only. Let a G (0, 1) be flxed in the remainder of this work. 
We define the set of admissible functions to be 

V := {/ e /i2+2-gM7)+-(g) ; |l/|l^(s) < 1}, 

where § is the unit circle. Note that when we neglect surface tension V C /i^"^"(S), 
whereas for 7 > we have V C /i^^"(§). The small Holder space /i™'+^(§), m € N 
and /3 € (0, 1), is defined as the completion of C°°(S) in C"^~^^{S), the classical Holder 
spaces. Furthermore, we set F^ := S x {k} for k € { — 1, 0, 1}. 

To incorporate time, we fix T > 0. If / G C{[0,T],V)r\C^{[0,T],h^+°'(S)) describes 
the evolution of the interface separating fluids, then, at any time t G [0,T], we have 
that ^±{t) = ^±(/(t)), where, given /i G V, we set 

n.ih):={{x,y)en : -1< y < /i(.x)}, 

n+{h) := {(x,y) G : h{x) < y < 1}. 

Consequently, it holds that F(i) = F(/(i)) for aU t G [0,r], with 

F(/i) := {(x,/i(x)) : X gS}, 

for all h (^V. Let us describe the kinematic condition (2.1e) in this particular context. 
Consider the evolution of a particle {x{t),y{t)) on the interface F(/(t)). Since particles 
on the interface separating the fluids remain there as time evolves, we obtain from 
F(i) = F(/(i)), that y{t) = f{t,x{t)) for all t G [0,r]. Differentiating this equation 
with respect to the time variable t yields dtf = { — f',l){x',y'), hence, the normal 
velocity V satisfies the relation 
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We denote herein by /' the spatial derivative of /. Thereby, the local study of problem 
(2.1) reduces to the following system 

Au± = in n±{f{t)), 0<t<T, 

dyU+ = 02 on Fi, < t < T, 

U- = gi on r_i, < t < T, 

(2.2) 
u+-u^ = ^uru) + ^f o^ r(/(i)), < f < T, 



9i/ + A;/x±iyrT7^9.n± = on r(/(i)), 0<f<T, 

/(O) = /o, 

where /o € V determines the initial shape of Fq, and to simplify notation, we set 

w.= g{p+- p^). (2.3) 

Vice versa, defining v± by (2.1a), system (2.1) can be deduced from (2.2). The 
functions gi and g2 are assumed to satisfy 

gi G C([0, oo), /i2+"(S)) and g2 € C([0, oo), /ii+°(S)). 

A triple (/, u+,n_) is called a classical Holder solution of (2.2) on [0, T] if 

/eC([o,r],v)nci([o,T],/ii+"(S)), 
n±(i,.)e&yc2+-(j7±(/(t))), te[o,r], 

and if (/, n+,ti_) satisfies (2.2) pointwise. Given / G V, we let buc ^"(f^±(/)) de- 
note the completion of the smooth functions BUC^{ft±{f)) in the Banach space 
BUC'^~^°^{VL±{f)). The first main result of this paper is the following theorem: 

Theorem 2.1 (Local well-posedness). Let 7 € [0, 00), ci,C2 G M, and assume that 

7>0 or ff(p+-p_) + C2f^-lj <0. (2.4) 

Then there exist open neighbourhoods of the zero function Oi C /i'^°(S), i € {1,2}, 
and O CV, such that for all /q € O and gi G C([0, 00), q + Oi), i = 1,2, there exists 
^(/o) > and a unique maximal Holder solution (/ =: /(•; /o),ti+,ti_) of problem 
(2.2) on [0,r(/o)) which fulfills f{t) G O /or all t G [0,r(/o)). 
-(/ 9i > 52 «s 0/ class C^ , then 

{{t,fo) : /oGOandtG(0,r(/o))}^/i2+2-g-(7)+-(g)^ ^^j^-^^f^^.j^^ 

is of class C^ . 

Furthermore, i/7 = 0, gi,g2 o-fe constants and if 

g{p+-p.) + g2(^-l\>Q, (2.5) 

then the linearised problem is ill-posed in the sense of Hadamard. 

Remark 2.2. If 7 > then we may choose Oi = h^~^°'{S), i G {1,2}. This means 
that we need no restriction on the boundary data when surface tension effects are 
considered. 
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Notice that if gi and g2 depend only upon time and are constant in the spatial 
variable, then, for T > small enough. 



fit) :-- 



-— f 92is)ds, te[0,T], (2.6) 

^J'+ Jo 

is a solution of the problem (see also [23]). This is a flat interface which moves 
according to the sign of g2- Consider now the tuple {f,u^,U-) to be an arbitrary 
solution of (2.2). Stokes' theorem yields then 

-y-r: f dx = - -s/l + f^dyU^ dx = -27r / dyUj^ = -Itx / g^ dx, 

which shows that the upper fluid moves towards the bottom if g2 > 0. If 52 < 
the fluid on the bottom must rise. This behaviour may be seen very well also from 
relation (2.6). Moreover, if g2 has integral mean equal to zero, then the moving 
interface oscillates around some constant value, since the formula above yields in this 
case 

f dx = const. 

Summarising, the direction of the flow is determined only by the boundary data g2, 
while gi is related only with the values of the potentials u± . 

Remark 2.3. The result stated in Theorem 2.1 is consistent with previous results on 
this problem. For example, if the densities of the fluids are equal, but the more 
viscous fluid drives forward the less viscous one the problem is well-posed, cf. [21-23]. 
Furthermore, if the lower-viscosity fluid expands into the higher-viscosity fluid, then 
the linearised problem is ill-posed, a result which is related with the investigation 
in [21]. 

Our analysis discloses that the problem is well-posed also when />_ < />+ and 
^_ > /i+. In this case 

gfi+jp- - p+) 

C2 < < 0, 

thus the less dense fluid replaces the less viscous one, but the velocity of the flow 
must be sufficiently large. Moreover, when p_ > pj^ and ^_ < /i_|_, the problem is 
well-posed provided 

gp+{p- - p+) 

C2 > . 

P- - p+ 

Notice that in this situation C2 may be negative, which shows that the less viscous fluid 
may drive upwards the less dense one if the velocity of the system is small enough. It 
is known that water is more dense than oil but oil is more viscous so that it adheres 
more strongly to the pores of the medium, cf. [4]. Hence, this relation determines also 
an optimal value of the normal velocity at which the water may replace the oil in an 
oil reservoir when surface tension effects are neglected. Above this threshold value 
fingering may occur and water fingers could penetrate into the oil domain. In view of 
f2.1c) we obtain: 



ON THE PARABOLICITY OF THE MUSKAT PROBLEM 7 

Corollary 2.4. The optimal normal velocity at which oil is driven upwards by water in 
an oil reservoir, when surface tension effects are neglected, is 

K.. = ^'^"""""\ (2.7) 

/i+ - ^_ 

if we identify the fluid on the bottom to be water and the one above as oil. 

Relation (2.4) was found also before, of. [20], when studying the Saffman- Taylor 
linearised stability of the flat interface between two fluids in motion under small 
perturbations. 



3. The operator equation 

In order to prove Theorem 2.1 we first transform the system (2.2) into a problem 
on a fixed domain. We then use certain solution operators for elliptic boundary value 
problems associated to (2.2) to rewrite the original problem as an abstract evolution 
equation. Let us introduce some notation. We set Q± := $^±(0) and identify the 
boundary Fq = S x {0} with the unit circle S. For every / E V, we define the 
mappings (pf± = {(l)j^,(pj^) : fi± -^ ^±{f) by 



(t>f±(.x, y) := {x, 2/ + (1 T v)f{x)), {x, y) £ n 



±- 



Note that (j)f± are diffeomorphisms mapping the reference domains Q± onto il±(/), 
i.e. ^f± G Z?z#^+"(fi±,0±(/)). Moreover, </./±(r±i) = r±i, ^f±{To) = T{f), and the 
inverse mappings ipf^ are given by the expressions 

'>Pf±{x,y) ■.= <j)jl{x,y) = ix,———] , {x,y) G 0±(/). 

These diffeomorphisms induce push-forward and pull-back operators 

(^{^ : buc^+''{n±) -^ &uc2+"(Jl±(/)), v^vo V'/i, 



7± 



: 6wc2+"(J7±(/)) -^ 6iic2+°(fi±), u^uo </.y±, 



which are, in virtue of the mean value theorem, isomorphisms, i.e. {4>*f^) ^ = </>{ . 
We now use these operators to transform system (2.2) on our fixed reference domains 
Q±. The price to be paid for doing that is the fact that the differential operators 
we have to study are more involved. Given / G V, we define the operators A±{f) : 
buc'^~^°'{i}±) —T' buc°'{Q,±) by the relation 

^±(/):=</'/±A</.{±. 

The uniformly elliptic operators A±{f) depend analytically on the a priori unknown 
function /, i.e. 

A± G C7'^(V,/:(6uc2+"(0±), 6uc"(0±))), (3.1) 
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since we have the fohowing representations of A±: 

^ {lTf)f"±2f'^ d 
-l=Fy .. ^ .so 7^ for/GV. 

(1 T /)^ 9y 

Indeed, A±{f) is a hnear combination of differential operators of order less or equal 
to 2 with coefficients depending analytically on /. Furthermore, A±{f) have the 
following geometric interpretation. Let rj be the standard metric on M^. The diffeo- 
morphisms (/)^-|- induce Riemannian metrics, 4'*fj_T], on 0±, i.e. 

</'/±??L(e,C) := ^|<^^^(,)W/±(x)^,a0/±(x)C) = ei<l>f±ix)f dcPf ±ix)C 

for X G Q± and tangent vectors ^, C ^ Tr^^±. With this notation, A±{f) are exactly the 
Laplace-Beltrami operators corresponding to the Riemannian manifolds {^±,4>*f-i^r]), 
respectively. Furthermore, given / E V, we define the boundary operators B±{f) : 

B±{f)v± := fc^iHr,/.}±(V((/){±t-±)|(-/', 1)), v± G 6uc2+"(f]±), 

where tr is the trace operator with respect to S, i.e. trt'-i-(x) := f±(x,0) for all x G S 
and v± G BUC{i}±). A simple computation shows 

B±{f)v± = — -tr — / tr 



/i± V 1 T / dy dx 

so that Bj- depend analytically on /, i.e. 

B± G C"^(V,/:(6uc2+"(rj±),/ti+" (§))). (3.2) 

Finally, let B{f) : 6uc2+"([7+) -^ /ii+"(S) be given by 

B{f)v+ := ——:tTidyV+, v+ G 6tic2+°($7+), 

with tri the trace operator with respect to Fi, that is trif+(x) := ^4.(3;, 1) for all 
j; G S and Vjf- G 5{/C(r2+). Clearly, we also have 

B G C'^(V,/:(6uc2+"(0+),/ii+°(S))). (3.3) 

Letting v± = 4'*f±u±, we see that problem (2.2) is equivalent to the following system 

■A±{f)v± = in n±, 0<t<T, 

BU)v+ = 92 on Fi, 0<t<r, 

V- = q^ on F_i , < t < T, 

(3.4) 
u_l_ — u_ = 7k(/) + wf on Fq, < t < T, 

dtf + B±{f)v± = on Fo, 0<t<r, 

/(O) = /o, 
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where k : /i^+"(S) -^ h'^+'^{S), with 

is the pulled-back curvature operator. Obviously, k depends analytically on / as well, 
i.e. K G C"^(/i'^+"(S),/i2+"(S)). The notion of classical solution for (3.4) is defined 
similarly to that for (2.2). The systems (2.2) and (3.4) are equivalent in the following 
sense: 

Lemma 3.1. Let {f,u^,U-) be solution of (2.2). The tuple {f,(j)*r,u^,(p*r_U-) is then 
a solution to (3.4). Vice versa, if {f,v^,V-) is a classical solution of (3.4), then 
{f,4>i V-\-,cl)l~V-) is a classical solution of (2.2). 

Proof. The main difficulty lies in showing that, if / G V, then u± S buc'^~^°'(Q±(f)) 
exactly when v± G buc'^~^°' {^±) . However, this equivalence can be proved by using 
the mean value theorem as in [10, Lemma 1.2]. D D 

We emphasize that the mapping /, describing the moving boundary separating the 
fiuids, is preserved by this transformation. Although the transformation we made 
above has the draw-back of introducing additional nonlinear coefficients, it allows us 
to rewrite the problem as a nonlinear operator equation. 

3.1. The operator equation. We now define solution operators to mixed boundary 
value problems which are closely related to our system (3.4). Composing these opera- 
tors, we obtain thereafter an operator equation which is equivalent to the formulation 
(3.4). Moreover, the operators defined below will reveal that it is crucial to determine 
the mapping /, describing the evolution of the interface. The potentials v± are then 
the image of / under these operators. For this reason, we refer also only to / as being 
a solution of (3.4). 

Given / G V and {q,h) G /ii+"(S) x /i2+"(S), let T{f,q,h) G 6uc2+"(J7_) denote 
the solution of the linear, elliptic mixed boundary value problem 

A-{f)v- = in Vt_, 

B.{f)v. = q on To, (3.5) 

V- = h on r_i. 

Further on, we define 5 : V x h^+°'{S) x /i2+"(S) -^ buc'^+"'{n+) by writing S{f,q,h) 
for the unique solution of the problem 

A+{f)v+ = in n+, 

B{f)v+ = q on Fi, (3.6) 

v+ = h + -fK{f) + wf on Fq. 

With this notation we observe that (/, w_|_,f_) is a solution of (3.4) if and only if 
/(O) = /o, V- = T{f,-dtf,gi), v+ = S{f,g2,tT:v_), and 

dtf + B+if)Sif,g2,trTif,-dtf,gi)) = 0. (3.7) 
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We have obtained in this way an operator equation only for the mapping / which has 
to be solved. A possible approach to (3.7) is to use the implicit function theorem and 
Newton's iteration method, similarly as in [22]. We proceed differently by decompos- 
ing the operators T and S appropriately so that we may rewrite the equation (3.7) as 
an abstract Cauchy problem which has order 1 + 2sign (7). 

Let us start by noticing that T(/, q, h) = 7i(/)/i + 72(/)'Z, where, given / G V, the 
linear operators 7i(/) and 72(/) are defined by 

ri(/)/i:=(^_(/),^_(/),tr)-i(0,0,/i), r2(/)g:=(^_(/),S_(/),tr)-i(0,g,0). 

Similarly, we decompose 5(/, g, h) = Si{f)h + S2{f)q + S3{f), with 

Siif)h := iA+if),Bif),tr)~\0,0,h), S2{f)q := (^+(/), e(/), tr)-i(0,g,0), 

Ssif) := (A+if), S(/), tr)-i (0, 0, 7ac(/) + wf). 

In virtue of (3.1)-(3.3), and taking also into consideration that the operator mapping 
a bijective linear operator onto its inverse is analytical, we obtain that the linear 
operators %{/), 5j(/), i = 1,2, as well as 53(/), depend analytically on /. With this 
notation, problem (3.7) rewrites 

=dtf + B+{f)S{f, 52, tr TiU)gi - tr T2{mf)) 

= [idhi+.^s)-B+{f)Si{f)trT2{f)]dtf + B+{f)S3{f) + B+{f)Siif)tvTiif)gi 

+ B+if)S2if)g2. 

The new expression of (3.7) is more involved, but we remark that the time deriv- 
ative dtf appears in the equality above as the argument of an operator belonging to 
C{h^~^'^ (§:)) . This operator is invertible for all / G V. 

Lemma 3.2. Given / G V, the operator G{f) '■= id/ji+c«(§) — ;B+(/)5i(/)tr 72(/) is an 
isomorphism, i.e. G{f) € Isom{h}'^°' {E>)) . 

Proof. Let p,q ^ h^~^°'{S) be given such that G{f)q = P- If u^ '■= {T2{f)q) ° 4'J_ and 
u+ := (5i(/)tr 72(/)^) ° 4''f+j then {u-,u+) solves the following system 



(3.5 



Au± 


= 





in 


^±{f), 


n_ 


-- 





on 
on 


r-i, 


flZ^d,yU^ 


q/V^ + f' 


r(/), 


M_ — n+ 


= 





on 


r(/), 


dyU+ 


= 





on 


Ti. 


q = p is equivalent to 







kfi7.^d,u+ = {q-p)/y^TT7^ on r(/). (3.9) 
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Given p G C^'^°(S), consider the diffraction problem which is obtained by ehminating 
q from (3.9) and the fourth equation of (3.8) 

An± = in Q±{f), 

U- = on r_i, 

kfiZ^d,U--kfi^^d,u+ = p/y^TTJ^ on r(/), (3.10) 

u_ — ti+ = on r(/), 

dyU^ = on Fi. 

We infer from maximum principles for elliptic problems that, when p = 0, (3.10) has 
only the trivial solution {u^,u^) = (0,0). Whence, for general p G C^^"(S), (3.10) 
has a unique solution {u-,u+) G BUC'^+''{Q._) x 5{/C2+"(Q+) (see [16, Theorems 
1.61 and 16.2]). Consequently, we find from (3.9) a unique q G C^'''"(§) such that 
Q{f)q = p. Taking into consideration that Q depends analytically on /, the desired 
result follows from a density argument. D D 

Applying the inverse C/(/)~^ to the equation (3.7), we obtain an equivalent formu- 
lation of our original system (2.2): 

Proposition 3.3. A function f is a classical Holder solution of problem (2.2) if and 
only if f is a solution of the evolution equation 

dtf = ^t,f), /(0) = /o, (3.11) 

where $ : [0, oo) x V ^ h^+'^iS), <^{t, f) := $i(/) + $2(t, /), is defined by 

Mf)--=-{G{f)r'B+{f)Ss{f), (3.12) 

MtJ) ■■= - iGif)r' [B+{f)Si{f)tvTi{f)gi{t) + B+U)S2U)92{t)\ ■ (3.13) 

for all {t, f) G [0, cx)) x V. 

4. The proof of Theorem 2.1 

We use in this section the abstract theory for evolution equations as presented in [17] 
to prove the results stated in Theorem 2.1. Let us first emphasize that the operator $ 
is analytic with respect to / and the partial derivatives with respect to this variable are 
continuous in both variables (i, /). In order to prove the well-posedness of the problem 
(2.2), we are left to check that the Frechet derivative 9j$(0,0) generates a strongly 
continuous and analytic semigroup. The following lemma is, besides Proposition 3.3, 
the essential part in the proof of Theorem 2.1. 

Lemma 4.1. If gi{0) =: Cj G M, i = 1,2, then the derivative (9j$(0, 0) is a Fourier 
multiplier of order 1 + 2sign(7). More exactly, given f G /i2+2sign(7)+a(g) /g^ f = 
SmeZ o-mG-^^^ denote its Fourier series expansion. Then, we have 

9jcD(0,0)/= J;A„a„e^'"^ (4.1) 
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where 



fc|m| tanh(|?TT-|) 

Arr 



/i+ + /i_ tanh (m) 



C2 { l] +w - 7m^ 



meZ. (4.2) 



In the main body of this section we prove Lemma 4.1 and, at the end of the 
section, we combine our results to prove existence, uniqueness, and dependence of 
solutions to (2.2) on the initial data and time. To this scope, we first determine an 
expansion for the isomorphism ^(0). Given q € h^~^°'{'S>), we consider its Fourier series 
q = X^mGZ '^m^*™^- ^^ ^^^ [10]' ^^ obtain that T2{0)q has the following expansion 

r2(0), = -(! + , )ao + - E ( ^(,^^,-^.) - ^(,^^,-^) J"^^-"'-- 

meZ\{0} "^ '^ ' '' ^ ' > '' 

Furthermore, given h € /i^+"(S) with h = X^mg^ ^mC*"*^, it holds that 

,_^ / p-'^p'^y prrip-my \ 

'5i(o)/^ = E ;;^r-73^ + t^^ttt^ ^rne-\ (4.3) 



+ e-™ e™ + e- 



and so 



5,(0). 75(0), . ^ V (4:;^ . ,1;^) i5^„.„e 






imx 



We set by convention tanh(?n,)/?TT, = 1 if tti = 0. Summarising, we obtain that 

g{0)q =y (l + ttl. tanh\m)] amc'"''' (4.4) 






for ah Q' G /i^+"(§) with q' = XlmeZ "m^* 

We begin by studying the mapping ^i. Taking into consideration that 53(0) = 0, 
we obtain that 

9$i(0) = -{g{0)r'B+{0)dS3{0), (4.5) 

where for / G /j2+2sign(7)+a(g)^ ^^ ^^^^ dS3{0)[f] = (A,S(0),tr )-i(0, 0,7/" + w/). 
The last relation follows in virtue of .4+(0) = A and 9k(0)[/] = /". Reconsidering 
(4.3), we see due to (4.4) that 9<^i(0) is a Fourier multiplication operator with 

9M0)[f] = E i'^^^'^'i, lm|(^-7m^)a^e- (4.6) 

;^ /i+ + /i_ tanh (m) 

for / = XlmGZ '^'»^*"*^ ^ /j2+2sign(7)+Q^g^_ This representation is sufficient to obtain 
the well-posedness of the problem when considering surface tension effects. How- 
ever, when surface tension is neglected, or when we study the stability properties of 
equilibria, it is necessary to analyse more closely the partial derivative (9j$2(0,0). 
In the following we set gi{0) = Ci, i = 1,2, where Cj G M. By the chain rule we get 

dfMO, 0)[/] = - (e(0))-i B+{0)dS2{0)[f]c2 - {g{0)r' dB+{0)[f]S2{0)c2 
- d ((a(-))"') (0)[/]S+(0)52(0)c2 =: hf + hf + hf, 
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for / G /i2+2sign(7)+a(g)^ Tj^jg jg ^^^ ^ ^he fact that ri(/)ci = ci, 5i(/)ci = ci, 
relations which imply that the first term of $2(0, /) is the zero function: 

- {gU))"' B+{f)Si{f)Ti{f)ci = 0, V/ e V. 

We determine next expansions for the linear operators /j, i E {1, 2, 3}. To this scope 
we note that 

a^.(o)[/i = -2(1 ^ rf/'^ ± 2/^ - (1 „)/"|, 

38,(0)1/1 = ^ (±/tr I - /'tr A) , .„a 96(0)[/l = /tr.| 

for all/€/i2+2'^^sn(7)+«(s). 

Let / G /i2+2sign{7)+Q(§) i^g gi^gj^ g^j^j denote by / = EmGZ«me*"''' i**^ Fourier 
series. We begin determining the expansion corresponding to Iif. The function 
952(0) [/]c2 is the solution of the boundary value problem 

Aw; = -dA+{Q)[f]S2{Q)c2 in 0+, 
Oyu; = -9S(0)[/]52(0)C2 on Ti, 
u; = on Fq, 

where 52(0)c2 = C2y in [7+. We make the following Fourier series ansatz 

5cS2(0)[/]c2 = J]u;^(2/)e*'"^ 

meZ 

Plugging this expression into the system found above, we are left, after identifying 
the coefficients of e**"^', to solve the problems 

w'^ - nT?Wm = C2amtv?{y - 1) in < y < 1, 

w'^{l) = -C2am on Fi, (4.7) 

■Wm{0) = on Fq. 

The solution of (4.7) is given by 

Wm.{y) =C2amtanh(m) 020^ h C2am(l - y) 

for all m € N and y € [0, 1]. Whence, 

B+{<d)w = —duw = y — (mtanh(m) - l)c2a^e*™^ 

and we obtain, in view of (4.4), the following expansion for Iif 

J f sr fe(^tanh(m)-l) ^^^ 

ill = - } ^ 7 — 7^7 — 7C2a'me- 

^^fi+ + fi_ tanh (m) 

Since iS2(0)c2 = C2y in ri_|_, a simple computation yields 

^2/ = - E ^ ^2, , C2a^e'"^^ 

fTt /^+ + M- tanh (m) 
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Lastly, we look for the Fourier series of the function I^f. To simplify notation we 
set J- := Q~^ and put c := A;c2//x+. In view of J~{f)Q{f) = id/ji+Q(§), the chain rule 

yields dT{0)[f]c = - {G{0)y^ dg{0)[f] (^(0))"^ c, while (4.4) implies that g{0)c = c. 
From the definition of ^, we get by differentiation that 

dg{0)[f]c = - 9S+(0)[/]5i(0)trr2(0)c - B+{0)dSi{0)[f]tvr2{0)c 
- S+(0)5i(0)tr 5X2(0) [/]c =: Jif + J2/ + J3/. 

We start by analysing J3. Differentiating (3.5) with respect to /, we conclude that 
972(0) [/]c is the solution of the linear elliptic problem 

Aw = -9^_(0)[/]75(0)c = c/i_A;-V"(l + y) in 1^-, 

dyw = -fi^k-^dB^{0)[f]T2{0)c = fdy{T2{0)c) = cn-k-^f on Fq, 
w = on F_i, 

where we made use of the relation 72(0)c = C/U_A;~^(1 + y). Expanding 

w{x,y) = Y^ 'u;^(y)e*™^, {x,y) S Jl_, 

we get that Wm is the solution of the system 

'' w'^ — m'^Wm = —cfi-k'^amm^i^ + y) in — 1 < y < 0, 

w'^{0) = cfi^k-^Um, (4.8) 

w;^(-l) = 0. 

Since Wm{y) = cji^k^'^ami^ +2/)) we find that 

<9'72(0)[/]c|y=o = ^ c/i_A;~^a.;„e'' 



Jimx 



Together with (4.3), we then get 

,-^ / p-rripmy m -my \ 

5i(o)tr9r2(o)[/]c = Y. -i^^^-' {-^^ri—^^^^ri—^)-rne-^- 

Whence, we have shown that J3 is a Fourier multiplier of the following form 
J3 5^«me^"'^= Y ^^"^tanh(m)o^e^"^^ 

meZ meZ\{0} 

To find an expansion for J^f we note first that 72(0)c = /U-c(l + y\ which leads to 
5i(/)tr72(0)c = c[i- for all / G V. This implies that J2/ is the zero function. Since 
5i(0)tr 72(0)c = /-i-C, we finally get J\f = 0. We have thus shown that dg{0)[f]c = 
J^f for all / G /i2+2sign(7)+a(^g^ Summarising, we find the following relation 

N r „-, ^-^ /m- \ kmtanh(m) ,^^ 

dfMO,0)[f] = E — - 1 —- , J. , C2ame'^-, (4.9) 

;^ V/^+ J fi++fi- tanh (m) 

which leads, together with (4.6), to (4.1). 
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We are now prepared to prove the main result stated in Theorem 2.1. The ar- 
gumentation strongly relies on the regularity of $, relation (4.1), and well-known 
interpolation properties of the small Holder spaces 

(/i'""(S),/i'"i(S))e = /i(^~^)'"»+^'"i(S), (4.10) 

if G (0, 1) and (1 - e)ao + 9ai ^ N. 

Proof of Theorem 2.1. Assume first 7 = 0. Using [10, Theorem 3.4] we obtain that 
the Frechet derivative 9/^(0,0), which has order 1, generates a strongly continuous 
and analytic semigroup in C{h^~^^{S)) for all /3 € (0, 1), if gi{0) = Cj € M, i = 1,2, 
and (2.4) is satisfied. More precisely, with the notation used in [2], this is written as 

Pick some /3 < a. Since 'H{h'^~^^{S),h^~^^{'B)) is an open subset of the space of 
bounded and linear operators >C(/i^+^(S), /i^"'"'^(S)) we find by continuity open neigh- 
bourhoods of the zero function Oi C /i*'^"(S), i G {1)2}, and O C V such that 
-df^{t,f) € 'H(/i2+/5(S),/ii+/5(S)) for aU / G O and gi'{t) e c^ + d. The existence 
result follows now in virtue of [17, Theorem 8.4.1] and relation (4.10). The regularity 
assertion is obtained by using [17, Corollary 8.4.6]. 

When 7 > there are no restrictions on the smallness of ^fj, 1 < i < 2. This is 
due to the fact that the derivative df^{t,f) is a third order operator, and all terms 
not containing 7 are treated as lower order perturbations. In view of [17, Proposition 
2.4.1] we obtain the desired local well-posedness result. 

If 7 = 0, ^j G M, i = 1, 2, and (2.5) is fulfilled, then the linearised problem 

dth = df<^{0,0)h, /i(0)=0, 

is ill-posed in the sense of Hadamard in /i^"'""(S), since 9j$(0,0) does not generate a 
strongly continuous semigroup in this case. This completes the proof. D D 

5. Stability properties 

In this section we assume that gi = const, and g2 = 0. Then, the flat interface /* = 
is, in view of (2.6), a steady-state solution with the constant pressure distributions 
U- = ii+ = gi- Note that in this case we deal with an autonomous equation since 
the operator ^ in (3.11) does no longer depend on the time variable t ($2 =0). In 
the following we are interested in the stability properties of the equilibrium f^, . Before 
doing this, we observe that: 

Remark 5.1 (Conservation of volume). Let 5^1 G M and g2 = 0. If / is a solution of 
problem (2.2) then 



Js 



f{t,x) dx = const. 

as long as the solution to (2.2) exists. 

Proof. See the discussion below relation (2.6). D D 

Following Remark 5.1, when studying the stability properties of the equilibrium f^, 
we are led to the assumption /o G /iq ^^^"^'~'' "(S), where given m G N, the space 
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/i™ (S) consists of the functions in /i™"^"(S) having integral mean equal to 0. Hence, 
we have to restrict our problem to the set Vq := V PI /iq''~"(S). 

Lemma 5.2. Given f G Vq, we have that ^{fo) G /iq'^"(S). 

Proof. We begin by showing that Q{f) € >Cis(/iQ^"(S)) for all / € Vq, which is equiv- 
alent to showing that B+{f)Si{f)tiT2{f)q € /io+"(§) for ah q € /io^"(§) and / e Vq. 
Indeed, setting u^ := Si{f)ti72{f)q, we have 

B+{f)Si{f)tiT2{f)qdx= l\-f'ix),l)-S/{cptv+)ix,f{x))dx 



k 



1 

2^ 



<9v (</'{"^w+ ) ds 
) ^ ^ 

Lastly, we still have to verify that B+{f)S3{f) G /iq~^"(S) for all / G Vq. However, 
this last statement follows by repeating the integration steps presented above. The 
desired assertion is now obtained in virtue of (3.12). D D 

We have thus reduced our problem to an autonomous evolution equation 

dtf = $(/), /(O) = /o, (5.1) 

where $ : Vq C /io+^'^^°^^^+"(§) -^ /iJ+"(S). From (4.1) and (4.2), we see that the 
spectrum of Frechet derivative 5^(0) G C{hQ ^'§° '■'>''' "(§), /ij!j+"(S)) consists only of 
the eigenvalues 9<l>(0) = {Am : 1 < m G N}, where Am are given by (4.2) with C2 = 0. 
In virtue of [17, Theorem 9.1.2] and the equivalence of problems (2.2) and (5.1) for 

initial data /o G h™^"" ^^^+'' {§) we get: 

Theorem 5.3 (Exponential stability). Assume that 

g{p+-p-)<l. (5.2) 

Then, the flat equilibrium /* = is exponentially stable. More precisely, given u G 
(0, tanh(l)A:(7 — w)/(fj,^ + /i_ tanh(l))), there exist positive constants 5, M such that 
for all /o G /iq ''^siH7J "^g^ ^^^^^ ||/o||;j2+2sign(7)+Q(-g) < 6, the solution to (2.2) exists 
in the large and 

ll/(0lk2+2sign(7)+Q(§) + \\dt f {t)\\f^l + a (^^) < Mc ^ 1 1 /q 1 1 /j2 + 2sign (7)+(i (§) 

for all t > 0. 

Moreover, if g{p^ — p„) > 7 > 0, this stationary solution is unstable. 

Remark 5.4. Let us notice that if surface tension effects are considered, then the 
equilibrium /* is stable also when the heavier fluid lies above, provided that the density 
jump across the interface is small compared with the surface tension coefficient, cf. 
(5.2). If we neglect the surface tension, then the equilibrium is stable only if the 
heavier fluid occupies the lower region of the cell. Moreover, this is also the only 
possible case which can be analysed, since via (2.4) we must have > 5(p+ — p~) 
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to ensure local well-posedness in this case. In the unstable case g{p+ — p) > 7 > 
pattern formation is evidenced in [14] by means of numerical simulations. 

6. Steady-state fingering patterns and instability 

In this last section we have a closer look at the stationary solutions of problem 
(2.2) under the same constant boundary conditions, gi = const, and §2 = 0, as in 
the previous section. We still assume that the cell contains equal amounts of both 
fluids, meaning that we are interested in determining the stationary solutions of the 
autonomous problem (5.1). First of all, we notice that if / G Vq is a stationary 
solution of (2.2), then the potentials u_ and u+ are both constant. This is due to 
the fact that they are both solutions of elliptic problems with homogeneous Neumann 
boundary conditions. We are led by the fifth equation of system (3.4) to the problem 
of determining the functions solving 

/" f 

^Tn~FJ2vU2 +^f = const, and / f dx = 0. (6.1) 

Recall by (2.3) that w = g{p+ — P-) is the constant which measures the density jump 
across the interface separating the fluids and parametrised by the function /. If the 
density of the fluid on the bottom of the cell is greater or equal then that of the fluid 
above, then (6.1) has only the trivial solution /=„ = (see e.g. [11]). When ro > 0, we 
shall use a bifurcation argument with the surface tension coefficient 7 as bifurcation 
parameter and obtain infinitely many global bifurcation branches consisting only of 
stationary solutions of (2.2). Of course, in this situation we consider the operator 
4> in (5.1) to depend also on 7, i.e. $ = $(7,/). However it turns out to be more 
convenient to treat (6.1) instead of the operator equation $(7, /) = in /iq (S). 

We fix ro > and look for (7,/) € (0, 00) x C^(S) with ||/|jc{S) < 1 which solve 
(6.1). Therefore, we introduce an operator which enables us to consider both equations 
of (6.1) at once. Since all the solutions of (6.1) are smooth, for even functions, the 
problem of finding the solutions of (6.1) is equivalent to determining the solutions of 
the equation 

T(7,/) = 0, (6.2) 

in (0,00) X U, where U := {f e C^,t"(S) : ll/llc(§) < 1} and the operator T : 
{0,00) xU CRx Co+"(S) -^ C^ddi^) is by definition the derivative of the left hand 
side of the first equation of (6.1) 

^(T' •^) ■= ^ (l + //2)3/2 - 3^ (l + ;/2)5/2 + ^f 

for (7, /) G (0, 00) X U. The space Cq g"(S) is the subspace of C^+"(S) consisting only 
of even functions with integral mean 0, and analogously C^^^(§) consists only of the 
odd functions in C"(§). Clearly, T depends analytically on its variables and T(7, 0) = 
for all 7 G (0, 00). Its Frechet derivative 9/T(7, 0) is a Fourier multiplication operator 
with 



5/T(7,0) 



^ Qm COs{mx) 



,m=l 



y (7m — wm) am sm.{mx) (6.3) 



m=l 
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Figure 2. Steady oil finger penetrating water. 

for all / = J2m=i «m cos(mx) G C|^+°(S). Given / G N \ {0}, we set 

7; := ml-\ (6.4) 

The first result of this section is the following global bifurcation theorem, which states 
that a global bifurcation branch emerges from the trivial fiat solution {(7, 0) : 7 > 0} 
at (7;, 0) for all 1 < Z G N, where 7^ is defined by (6.4), provided zu > 0. 

Theorem 6.1 (Steady-state fingering solutions). Let w > and I > 1. The point (7;, 0) 
belongs to the closure TZ of the set ofnontrivial solutions of (6.2) in (0, 00) xU. Denote 
by Ci the connected component of TZ to which (7;,0) belongs. Then Ci is unbounded in 
(0,00) xU. 

Additionally, Ci has, in a small neighbourhood 0/(7^,0), an analytic parametrisation 
ilhfi) ■■ i-Si,6i) -^ (0,00) xU, and 



fi{e) = ecos{lx) + 0{e'^), for e 



0. 



Moreover, any other pair (7, 0), with 7 > 0, is not a bifurcation point. 

Possible stationary fingering pattern solutions of problem (2.2) are pictured in Fig- 
ure 2. That C;, / > 1, is unbounded in (0, 00) x U means that either Ci is unbounded 
in M X C^"'""(S), or that Ci reaches the boundary of (0, cxd) x lA. It should be mentioned 
that both situations may occur [9]. 

Remark 6.2. For all/ G N, / > 1, we have that 7j'(0) = and 7^(0) > 0. Consequently, 
the bifurcation is supercritical. The bifurcation diagram is pictured in Figure 3. 

Concerning the stability properties of these finger-shaped steady-state solution we 
state: 

Theorem 6.3 (Instability of the fingering patterns). Let 1 < / G N 6e given. The 
stationary solution fi^e) of problem (2.2), when 7 = 7i(e), is unstable provided ^ \e\ 
is small enough. 

Proof of Theorem 6.1. It is well-known that bifurcation may occur at (7, 0) only if the 
derivative 9/T(7, 0) is not an isomorphism, which in view of (6.3), leads to 7 = 7; for 



ON THE PARABOLICITY OF THE MUSKAT PROBLEM 



19 



k 






Figure 3. The bifurcation diagram. 



some 1 < / G N. One can easily verify the assumption of the theorem on bifurcations 
from simple eigenvalues due to Crandall and Rabinowitz [7], as stated in [3, Theorem 
8.3.1], and obtain that analytic bifurcation branches {'yi,fi) : (—5,5) — > (0, oo) x U 
consisting entirely of steady-states solutions of (2.2), emerge at (7^, 0) for all 1 < / G N. 
In order to study the global behaviour of these branches we rewrite equation (6.2) 
as follows 

r - T^ + -/'(I + rr'' = o, (7, /) e (o, 00) x u. 



+ -/'(I + f'f^' 

7 



l + /'2 

We use now the property of the operator A : Cq^"{§) -^ C^ddi^)^ ^f = /'"' ^° ^le 
an isomorphism. By applying the inverse of A to the equation above we obtain an 
equivalent formulation for the problem (6.2) 

F{j,f) = f + H{j,f)=0, (6.5) 

where H : {0,oo)xU C Mx Cq"^"(S) -^ Cq^"(S) is the completely continuous operator 

3/7"' 



H{^J) = A'- 



w 



+ -/'(I + ry 

7 



l + /'2 

Using Fourier expansions for functions in Cg'^"(S) as we did before, we conclude 
that the derivative dfF('y,0), has an odd crossing number at 7 = 7; for all 1 < I, 
cf. [18, Definition II. 3.1]. Particularly, the index i{dfF{'y,0),0) jumps at 7 = 7^ from 
1 to —1, or vice versa. We infer from the global Rabinowitz bifurcation theorem [18, 
Theorem II. 3. 3] that the connected component Ci of 7^ to which (7;,0) belongs, is 
either unbounded in (0, 00) x U or contains some other bifurcation point (7^, 0), with 
k ^ I. Using methods from the theory of ordinary differential equations, it is possible 
to show [9, Theorems 3.2 and 4.1] that the second alternative cannot occur, that is 
all branches Ci are unbounded in (0, 00) x U and pairwise disjoint. 

The derivative 7^(0) is zero due to the fact that T(7, — /) = for all (7,/) € TZ. 
Particularly, we have that fi{—e) = —fi{£) and 7i(e) = 7/(— e) for all e G (0,5). We 
compute now the second derivative 7;"(0). To this scope, for fixed 7 > 0, we define 
: M'^ — ^ M by the relation 

,2 



{x,y,z) =7 



(l + x2)3/2 



37- 



xy- 



l + x2)5/2 



+ mx for (x, y, z) G 
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It holds that T(7, /) = (j){f', /", /'") for ah f eU. Since dfjCpiO) = for ah 1 < i, j < 3 
we get that 9?rT(7,0) = 0. Furthermore, the third order partial derivatives of in 
all vanish, except for 

dfi3(l>{0) = -37, df^^^O) = -67, 

so that 9|^^T (7,0) [/,/,/] = -97/'V'" - I87/7'" for ah 7 > and / € C3+"(S). 
We may write for the analytic parametrisation //(e) = ecos(/2;) + T;(e), where ti{0) = 
r/(0) = 0, and T;(e) belongs to the closed complement of Ker5jT(7;,0) in Cq"^"(§). 
Differentiating the relation T(7;(e),ecos(/x) + T;(e)) = three times with respect to 
e, at e = 0, yields, in view of 7;(0) = and 5? rT(7;, 0) = 0, that 

9|;^T(7z, O)[cos{lx)f + 37r(0)52^T(7„ 0)[cos(Zx)] + 9jT(7„ 0)[r'"(0)] = 0. 

However, dfT(^i,0) has values in the complement of M • sin(/x) in C"(S), so that, by 
multiplying the relation above by sin(/x), followed by integration over the unit circle, 
we get 

^ l (g|^^T(7;,0)[cos(/x)]3|sin/x) ^ 37,^ 
^'^ ^ 3 (a2^T(7„0)[cos(/x)]|sin/x) 4 > ' 

where {f\g) := J^fgdx, f,g& L2(S), is the scalar product on -L2(S). In view of 7; = 
w/P, we conclude that 7^(0) = Sw/A for alH > 1. This finishes the proof. D D 

Finally, we come to the proof of Theorem 6.3. As in the previous section, the 
constant boundary conditions on F-i-i ensure that the operator $, defined by (3.11) 
(and which depends now also on 7), is independent of time t. The existence of the 
local bifurcation branches {■~fi,fi) can be also obtained by applying the theorem on 
bifurcations from simple eigenvalues due to Crandall and Rabinowitz to the operator 
equation 

$(7,/) = 0, (7,/)€(0,oo)xVo,e (6.6) 

where letting /i™+"(§) = {/ S /i'"+"(S) : / is even } n /i™+"(S), we set Vo,e := 
Vn /igg (S). A maximum principle argument shows that <I>(7,/) is even if / G Vo,e; 
meaning, in virtue of Lemma 5.2, that <I> : (0, 00) x Vo,e -^ ^oe (^) i^ ^ well-defined 
mapping. However, it is much more difficult to determine the derivative 7^(0) when 
applying bifurcation theory to the operator equation (6.6). 

When proving the instability of the steady-state solution (7i(e), //(e)), £ € {—5i,6i), 
1 < I, for the problem (5.1) it is more accessible to show that they are unstable steady- 
state solutions of the restriction 

dtf = H7,f), /(0) = /o, (6.7) 

where $ : (0, 00) x Vo,e — ^ ^o'e"(^)- Therefore we must study how the spectrum of the 
linearised operator varies along the bifurcation curves. Our main tool is the exchange 
of stability theorem of Crandall and Rabinowitz [8]. According to [8], we introduce 
first some notation. 
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Definition 6.4. Let E, F be Banach spaces and T,K ^ C{K,¥). Then /i € M is a 
K— simple eigenvalue ofT if 

diniKer (T — fiK) = codim (T — fiK) = 1, 

and, if Ker (T — i^lK) = span {u*}, then Ku^ ^ Im (T — fiK). 

The result proved by Crandall and Rabinowitz in [8] reads as follows: 

Theorem 6.5 (Crandall-Rabinowitz). Let E,F he Banach spaces, K G >C(E,F) and 
J- : (0, oo) xE ^^ F of class C^ with J^(7, 0) = near^^,. IfO is a K— simple eigenvalue 
ofT:= 5jJ-'(7*, 0) and a 9^rJ-'(7^,, 0)— simple eigenvalue ofT, then there exists a local 
curve (7(e), /(e)) such that (7(0), /(O)) = (7*,0) and F{-~^{e), f{e)) = Moreover, if 
-^(7>/) = 0, / / 0, and (7,/) is close to (7*,0), then {j,f) = (7(e), /(e)) for some 
e/0. 

Furthermore, there are real numbers /i(e), A(7) and vectors n(e), 77(7) S E suc/i i/iai 

9/ J-(7(e), f{e))uie) = iJi{e)Ku{e), dfFi^, 0)v{^) = X{^)Kv{j), 

with fj.{0) = A(0) = and u{0) = v{0). Each curve is C^ , with 

A'(7.)7^0 and lim -'^'^f'^^*^ = l. 

Proof of Theorem 6.3. In virtue of (4.1) and (4.2) we have 

f, / s A;777-tanh(m) , 9, 

(^ /i+ + /i_ tanh (m) 

Since Ai(7;) > for all / > 2 we conclude that the spectrum of the linearisation 
df^{'jl{e), fi{e)) contains positive eigenvalues if e is small, and the case Z > 2 is 
proved. 

If / = 1, then we find ourselves in the critical case when is the only non- 
negative point in the spectrum of 9/$(7i,0), since Ai(7x) = 0. The proof relies 
now strongly on Theorem 6.5. Recall that by (4.1)-(4.2) the derivative 0/-$(7, 0) € 
/:(/i^+-(S),/iJ+°(S)),isgivenby 

00 00 

5/$(7,0) ^ amCos{mx) = ^ Am(7)am cos(mx), 

m=l 771=1 

for all / = X^^=i am cos{mx) S /iq'|;"(S), with Am(7) as above. Setting T = df^{ji, 0) 
and I'C to be the inclusion /iq~^"(§) "^^ /iq~^°(§) we find all assumptions of Theorem 
6.5 fulfilled. Moreover, in virtue of X'{ji) < 0, and since 7i(e) > near e > (re- 
spectively negative when e < 0) (cf. Theorem 6.1), we conclude that the eigenvalue 
//(e) of 9j<I>(7i(e), /i(e)) must have positive sign for e close to 0. Hence, /i(e) is an 
unstable steady-state of (6.7) when 7 = 71(e). The proof is now completed. D D 
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7. Appendix 

Let us now assume that both fluids have the same viscosity fi = fj,±. We study the 
problem (2.2) in a frame moving with constant velocity 1/ > 0, by using the results 
already established in Theorem 2.1 and Theorem 5.3. More explicitly, in a infinitely 
long vertical column consisting of two fluids we describe the motion of the interface 
separating the fluids and which is located between two parallel lines r-|-(t) which 
move with constant velocity V. We impose constant normal velocity V on r_|_(t) and 
prescribe the pressure on r_(i). This situation is similar to that studied in [20,21] 
where it is assumed that the velocity at ±00 is equal to (0, V). 

Consider now the global solution / of problem (2.2) under constant boundary con- 
ditions (72 = and 5(1 = c € M, with initial data /o E /Iq~''"(S). We then define 

h{t) :=/(i) + ty e/i^+"(S); 

0+(t, h) := {(x, y) : h{t) <y<l + tV}; 

Q_(t, h) := {(x, y) : -1 + tV < y < h{t)}; 

I \ I TAN M^ , ^^ 

v±{t, X, y) := u±{x, y - tV) —y ± — -f , 

k I 

Setting r±(t) := S X {±1 + W}, we see that {h,Vj^,v. 
Muskat problem in the moving frame: 

Av± = 



--dyVj, = V 



0-1^ + °" 



{x,y) 


eQ±{t,h). 




is the global solution 


of the 


in 


n±{t,h), 




on 


r+(t), 




t on 
on 


r_(t), 
r{h{t)), 


(7.1) 



k 2 

v+ - V- = 7Kr(h,) + ^h 

d^h + ^^^—^d^v± = on T{h{t)), 

A* 

HO) = /o, 

with t G [0, 00) and C := {c+fiV)k. Defining the velocity field v^ by Darcy's law v^ := 
— {k/n)Vv± in Q±{t,h), the first boundary condition in (7.1) asserts that the normal 
velocity on r+(t) is constantly equal to V. Complementary to [21, Theorem 4.2] (we 
treat here fluids with different densities and the same viscosity), our next result states 
not only global existence for small data but also exponentially fast convergence to a 
flat solution. In virtue of Theorem 2.1, Theorem 5.3, and Theorem 6.1 we obtain: 

Theorem 7.1. Let (5.2) hold true andC G M, y > be fixed. Givenu € (0, tanh(l)(7 — 
tu)/(//(l + tanh(l)))), there exist positive constants 6,M such that for all fo € /iq'''"(S) 
with ||/o||/i4+Q(§) < S, the unique solution to (7.1) exists in the large and 

\\h{t) - iy|U4+.(§) + \\dth{t) - y|Ui+c.(§) < Me-'^*||/o|U4+.(§) 
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for allt > 0. Moreover, if p+ > p- andj > 0, there exist traveling finger shaped global 
solutions of (7.1). 

Remark 7.2. Though the normal velocity on r_(t) is not constant, in the limit we 
have {v_\i') -^t-^oo ^ exponentially fast if (5.2) holds true. If 7 = 0, relation (5.2) 
ensures that the more dense fluid lies beneath. 

Acknowledgement. The authors wish to thank the anonymous referees for their con- 
structive suggestions and comments which have improved the quality of the paper.. 
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